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Contributions of Spherical Harmonics 

to 

Magnetic and Gravitational Fields 

Carlos Roithmayr 


Introduction 

Gravitational forces are of cardinal importance in the dynamics of spacecraft; magnetic 
attractions sometime play a significant role also, as was the case with the Long Duration 
Exposure Facility, and as is now true for the first segment of Space Station Freedom. Both 
satellites depend on gravitational moment and a device known as a magnetic damper to 
stabilize their orientation. 

Magnetic fields are mathematically similar to gravitational fields in one important re- 
spect: each can be regarded as a gradient of a potential function that, in turn, can be 
described as an infinite series of spherical harmonics. Consequently, the two fields can be 
computed, in part, with quantities that need only be evaluated once, resulting in a savings 
of time when both fields are needed. 

The objective of this material is to present magnetic field and gravitational force ex- 
pressions, and point out the terms that belong to both — this is accomplished in Sections 1 
and 2. Section 3 contains the deductive reasoning with which one obtains the expressions 
of interest. Finally, examples in Section 4 show these equations can be used to reproduce 
others that arise in connection with special cases such as the magnetic field produced by a 
tilted dipole, and gravitational force exerted by an oblate spheroid. The mathematics are 
discussed in the context of terrestrial fields; however, by substituting appropriate constants, 
the results can be made applicable to fields belonging to other celestial bodies. 

The expressions presented here share the characteristics of algorithms set forth in 
Refs. [1]— [4] for computing gravitational force. In particular, computation is performed 
speedily by means of recursion formulae, and the expressions do not suffer from the short- 
coming of a singularity when evaluated at points that lie on the polar axis. 
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1 Contributions to Magnetic Field 


The magnetic field vector B at a point Q above the Earth’s surface can be expressed as the 
sum 

oo n 

B = E E B n,m (1) 

n= l m = 0 

where B n ,m represents the contribution to B of the spherical harmonic of degree n and 
order m, and is given by 


B 77 ,. m — 


Kn,m a n+2 f 9n ,m,Cm T h n m.Sr 

j^n+m+1 | 


R 


[(s A A„ im+ i + (n + m + 1 )A n , m ) r - A n , m+ ie 3 ] 


TflA n ,m \i,9n,mPm— 1 T "I" {,h n ,mPm— 1 9n,mSm— 1 )e 2 ]} 


(2) 


Here, a is the mean radius of the Earth (6371 km: p. 25 , Ref. [ 5 ]) adopted for the Interna- 
tional Geomagnetic Reference Field; R is the magnitude of R, the position vector from the 
geometric center of the Earth to Q; r is a unit vector in the direction of R; and g n ,m and 
h n ,m represent Gauss coefficients of degree n and order m. Mutually perpendicular unit 
vectors ei, 62, and £3 are fixed in the Earth: ei lies in the equatorial plane parallel to a line 
intersecting Earth’s geometric center and the Greenwich meridian, £3 is in the direction of 
the north polar axis, and e2 = £3 x ei. 

The coefficients K n m in Eq. ( 2 ) are given by 

K n ,o = 1 (n = 1 , . . . , 00) ( 3 ) 

whenever m = 0 (in addition, K 1.1 = 1 ); otherwise, they can be obtained with one of two 
recursion formulae: 

YY}, \ 2 

— - — - K n _i m [m= l,...,oo; n> (m + 1)] (4) 

n + mj 

_ 1 

K n ,m = [(n + m)(n - m + 1)] 2 K n , m -\ (m = 2, . . . , 00 ; n > m) (5) 

A n m and A n , m+ 1 are referred to as derived Legendre polynomials of degree n, and orders 
m and m +1 respectively. The notation A njm (u) is used to indicate that A n ^ m is, in general, 
a function of an argument u ; however, for the sake of clarity, the arguments of A n ^ n and 
A„,m+ 1 do not appear explicitly in Eq. ( 2 ). When n = m , the derived Legendre polynomial 
is given by a simple expression that can be written in recursive form. 

A n , n = ( 1 )( 3 )( 5 ) • • • (2n - 1) (n = 1, ... ,00) 

= ( 2 n - l)A n _i, n _i (n = 2 ,..., 00) (6) 

It should be noted that A2 2 j-- - are not functions of an argument; therefore, they 

need only be computed once to provide the seeds for a recursion relation shown in Ref. [6] 


2 



to give numerically stable and precise results: 

A n , m (u) = — - — [(2n - l)uA n _ l m (u) - ( n + m - l)A n _ 2m (u)\ 
n — m 

[m = 0, . . . , oo; n > (m + 1)] (7) 

A ntm (u) vanishes when m > n, and Ao } q(u) = 1. Arguments of the derived Legendre 
polynomials in Eq. (2) are sin A (shortened to sa), where A is the geographic latitude of Q. 
s A = r • e 3 . 

The symbols S m and C m appearing in Eq. (2) will be defined presently. Their values for 
m = 0 and m = 1, 

So = 0, C 0 = 1 (8) 

Si = R • e2, Ci = R • ei (9) 

provide a means of starting the recursion relations given on p. 4 of Ref. [2] . 

Sm = SiC m _i -f- CiS ro _i, C m = C\C m —i — SiS m _i (10) 


2 Contributions to Gravitational Field 


Like the magnetic field vector, the gravitational force per unit mass f exerted by the Earth 
at a point Q above the Earth’s surface, can be regarded as a sum: 


f = 



+ EE* 

n=2 m = 0 


(ii) 


where /a, the gravitational parameter, is the product of G, the universal gravitational con- 
stant, and M, the Earth’s mass; the position vector from Earth’s mass center to Q has 
magnitude R and the direction of unit vector r ; and f n m represents the contribution to f 
of the spherical harmonic of degree n and order m, given by 


f — 

L n.. m. — 


f C nt mC m T S n ,m<Sm 


ftn+m+1 


R 


[A n>m+ ie 3 - (saAj ;TO+ i + (n + m + 1 )A n , m ) r] 


+ mA n , m [(C n ^ m C m— 1 Sn,mSrri— 1 )^1 (Sn,mCm— 1 C n jnS m — 1 )^ 2 ] j 1 (12) 


R® is the mean equatorial radius of the Earth (6,378,139 m), and C n ^ m and S U)m are 
gravitational coefficients of degree n and order m. Mutually perpendicular unit vectors 
ei, e 2 , and e 3 have the same directions as in Eq. (2). A nm , S in , and C m have the same 
meanings as in Eq. (2); consequently, they need only be computed once when constructing 
magnetic and gravitational fields. 

It is not too difficult to show that Eq. (11) and Eq. (12), which has the form f n , m = 
x n , m & l + yn,m & 2 + 3 + w n .r r A [Eq. (2) has this form as well], are in concert with 

Eq. (31) of Ref. [1], whose form is f = Xei + Te 2 + Ze^ + Wv: in other words, W = 
-y 'JR 2 + E“=2 Em=o w n,m, X = E^= 2 Em=o x n,m, and so forth. (The sum to the right of 
the second = sign in Eqs. (30a) of Ref. [1] should be preceded by a minus sign.) 
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3 Derivations 

Eq. (2) is similar in form to Eq. (12) because B and f are each obtained by differentiating 
a potential function that is expressed in terms of spherical harmonics. A potential function 
associated with the magnetic field is given below, followed by steps in the differentiation 
process that results in Eq. (2). Next, a gravitational potential function is displayed. The 
differentiation process that leads to Eq. (12) is virtually identical to that used with the 
magnetic potential; therefore, it does not bear repeating. 


3.1 Magnetic Field 

The magnetic field vector B is related to Vp, a scalar geomagnetic potential, by the equation 

B = -VV/3 (13) 


where V denotes differentiation with respect to R: B is the negative of the spatial gradient 
of Vp. Section 2.9 of Ref. [7] contains a thorough explanation of how one differentiates 
with respect to a vector. The potential Vp satisfies Laplace’s equation and, consequently, 
is expressed in Eq. (2.13) of Ref. [5], and in Eq. (9.0,6) of Ref. [8], as an infinite series of 
spherical harmonics. 


Vp = a ^ 


Au) 


n+1 


53 P n( cos 8) 9n,m cos m(f) + h n;m sin mcj) 


m = 0 


(14) 


where a, R, g n ,m > and h n , m have the same meanings as in Eq. (2), and <fi is the geographic 
longitude of Q measured Eastward from the Greenwich meridian. P™(cos0) are Schmidt 
functions (they will be discussed directly) whose arguments are cos 9, where 9 is the geo- 
graphic co-latitude of Q. Now, Eq. (14) can be rewritten as 


Vp = E E V n,m (15) 

n= 1 m = 0 


so long as V n{m is defined as 


A 


a \ 


n+1 


V n ,m = a ( ^ J cos 9) \g n ,m cos mcj) + hn y m sin mcj) 


(16) 


For the sake of clarity, we dispense with the subscript (3 for V n , m . Thus, in view of Eqs. (13) 
and (15), B can be expressed as in Eq. (1) 


b = —Wp = -v 53 53 v n , m = 53 53 B n . m (u) 

n=l m=0 n=l m=0 


with the contribution to B of the spherical harmonic of degree n and order m defined simply 
as 


B, 


= -VK 


(18) 
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The Schmidt functions P™ are of degree n and order m, and are related [Ref. [9] , Chapter 
9, Eq. (13)] to P n .m, associated Legendre functions of the first kind. 


P™ (cos <9) 


0 (n— m)! l 2 
(ra+m)! 

P n ,m(cOS0) 


Pn,m(cos0) 


( m = 1, . . . , n) 
(m = 0) 


(19) 


Angle 0, the co-latitude of Q, is the complement of A, the latitude of Q ; therefore, the 
argument cos0 of P™ and P nm can be replaced with sa, which is shorthand for sin A. 
Products of P n ,m( s a) with cos mcj), or with sinm<^>, are referred to as tesseral harmonics of 
degree n and order m: they are functions of latitude and longitude. Tesseral harmonics of 
order zero are called zonal harmonics, and are independent of longitude. 

In the interest of economy in written symbols, the coefficients of P re , m in Eqs. (19) can 
be labeled K nm [see Eq. (2)] and defined as 


K — 


c\(n— m)\ 2 
(n+m)! 

1 


(m = 1, . . . , n) 
(m = 0) 


(20) 


The recursive relationships in Eqs. (4) and (5) follow from the first of these definitions; 
Eqs. (3) from the second. 

The associated Legendre functions P n ,m are, in turn, related to derived Legendre poly- 
nomials A n _ m by Eq. (9b) of Ref. [6] (see also p. 3 of Ref. [2]). 

Pn,m(s\) = (COS \) m A n , m (s\) = A n , m ( S A ) (21) 

where p is related to R by the definition 


5 2 = (R 


ei) 2 + (R • e 2 ) 


(22) 


In words, p is the magnitude of the projection of R onto Earth’s equatorial plane. 
The symbols S m and C m appearing in Eq. (2) are defined on p. 4 of Ref. [2], 


A 


S m = p m sin m<t>, 


A 


Cm = p m cos mcj) 


(23) 


where cj) has the same meaning as in Eq. (14). 

Differentiation of I+ n with respect to R is facilitated by substitution from Eqs. (19), 
(20), (21), and (23) into (16). 


T/ - / a \ n+1 

Vn.m, — & I I -K-i 


j -~n,m ) ^.n,m ($a) 1 9n,m COS 7Ti(j) + h n ^ m sin 171(f) 


= K, 


n,m 


77+2 a 9n,mCm ^n,m^m 
^ n,m ^Ti+m+l 


(24) 


K n mi a i Sn,m , and Ki.m are independent of R, so obtaining an expression for the gradient 
VK,m becomes a matter of producing the four gradients V R~^ n+m+l \ VC m , and 
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Eqs. (2.9.7)-(2.9.9) of Ref. [7] provide the machinery with which all the required gradi- 
ents can be formed: 

VR = U (25) 

Vi? = r (26) 

Vf = i? _1 (U-rr) (27) 

where U is the unit dyadic. (When U is pre or post multiplied by any vector v, the result 
is v, that is, U • v = v • U = v.) The first gradient is fashioned easily with the help of 
Eq. (26). 

VR -(n+m+ 1 ) = _( n + m + 1 ) ir (n+m+ 2 ) f . (28) 

Referring to Eq. (10) of Ref. [6], one can write VA njm (s A ) = A n)m+ iVs A , express s A in terms 
of r and 63 (63 is independent of R) in order to expedite differentiation with respect to R 
by means of Eq. (27), 

Vs A = V(r • e 3 ) = i?~ 1 (U - rr) • e 3 = [e 3 - (r • e 3 )r]/i? = (e 3 - s A r)/i? (29) 

and obtain the second gradient. 

V A nm (s A ) = (<53 s A f)^ n , m+ i/i? (30) 

Finally, the remaining two gradients can be produced from Eqs. (23) with the aid of Eqs. (22) 
and (25), 

V S m = fn{S m — ibi + C m _ie2), V C m = m(C m — iei — iS m _ie2) (31) 

as is suggested by Eqs. (7) in Ref. [2], 

The foregoing material puts one into position to express B n ,m as in Eq. (2) by substi- 
tuting from Eq. (24) into (18), and then employing the right hand members of Eqs. (28), 
(30), and (31). 

3.2 Gravitational Field 

The force f in a gravitational field is analagous to the vector B in a magnetic field: both 
are spatial gradients of scalar potential functions. At a point Q above Earth’s surface, a 
particle of unit mass is subject to a gravitational force f given by 

f = W 7 (32) 

where is a scalar gravitational potential, and V denotes differentiation with respect 
to R, the position vector from Earth’s mass center to Q. Like Vp, V-,, satisfies Laplace’s 
equation and can be expressed as an infinite sum of spherical harmonics. Eq. (2.13.12) of 
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Ref. [7] is one such example, written in the form adopted as a standard by the International 
Astronomical Union: 




R 


R \ n n 

1 + urn £ Pn,m (sin A) I Cn.m COS m 4 > + 5 n , m sin m(j) 

K ' m= 0 


n=2 
oo n 


+ EE ^n,m 

n=2 m= 0 


(33) 


where /j, i?0, C n . m , and 5 n)m have the same meanings as in Eq. (12), and Eqs. (21) and 
(23) can be engaged to define V n _ m : 


A d n / 1 C! n rn C m -\- S n ^ m S m 

V n ,m ~ H-n® ^n,m ^ n+m+1 

Thus, Eq. (11) is the result of substituting from Eq. (33) into (32). 


CXJ / 1 f i 

f = v i + E E w „, m = -Af + x; E f. 


n=2 m= 0 


n=2 m=0 


(34) 


(35) 


and the expression for f n m in Eq. (12) follows from differentiation of the right hand member 
of Eq. (34) [which is of the same form as Eq. (24)] in the much the same way as B n ,m results 
from differentiation of the right side of Eq. (24). /i. Rq, C' n . m , and S U)m are all independent 
of R. 


4 Examples 

Although Eqs. (2) and (12) are intended to be evaluated numerically, rather than analyt- 
ically, they can be used to produce recognizable expressions in a few important special 
cases. In this Section, Eqs. (2) and (12) are exercised through applications to a magnetic 
field associated with a tilted dipole, to the gravitational force exerted by a planet whose 
mass is distributed symmetrically about its polar axis, and to the gravitational attraction 
of an oblate spheroid — all of which have been dealt with elsewhere. 

4.1 Magnetic Field of a Tilted Dipole 

The geomagnetic field can be modeled as a tilted dipole by adding the contributions of 
spherical harmonics of degree 1, and orders 0 and 1. According to Ref. [5], approximately 
90% of the magnetic field at the Earth’s surface is described by this model. The derived 
Legendre polynomials identified with the dipole are [see Eqs. (6) and (7)] simply 

4i,o(sa) = sa, Aip(sa) = 1, 4i i 2(sa) = 0 (36) 

and Eq. (20) yields the coefficients 

Ki , o = K hl = 1 (37) 
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Substitution from Eqs. (36), (37), ( 8 ), and (9) into (2) leads, for n = 1, m = 0, and m= 1, 
to 



and then, after recognizing that s\ = r • 63 , to 

Bi, 0 + B u = 3 [3(r • M)r - M] (40) 


which is in agreement with Eq. (H-22) of Ref. [10] when M, the terrestrial dipole moment, 
is defined as M = g\ y \e\ + 1&2 + 51 , 0 ^ 3 . 


4.2 Gravitational Force Exerted by an Axisymmetric Body 

An expression for the gravitational force per unit mass, exerted by an axisymmetric body, 
is given in Ref. [11], and can be reproduced from Eqs. (11) and (12). The tesseral harmonics 
vanish from the gravitational potential as a consequence of axial symmetry, leaving only 
the zonal harmonics. Working only with terms of order 0, and substituting from Eq. ( 8 ) 
into ( 12 ), we get 

OO OO t~) 77, 

f = f n ,o = -jp* ~ I 1 5Z {^,ie 3 - [s A A n ,i + (n + l)A n>0 ] r} (41) 

n = 2 n= 2 


where the zonal harmonic coefficients, C H: 0 , have been renamed: C n , 0 = — J n . The derived 

Legendre polynomials are defined in Eq. (14) of Ref. [1] as A njTn (u) = d m P n (u) / du m , where 
P n (u) is a Legendre polynomial of degree n, with argument u. Denoting the first derivative 
(with respect to u) of P n (u ) as P' n (u) leads to 


. 00 o n 

f = - /' 5Z J nj ^2 { p n®3 - [s \P' n + (n + 1 )P n ] r} (42) 

n= 2 

which, upon use of the recursion relation P' l+1 (u) = uP' n {u) + (n + 1 )P n (u), becomes the 
expression given on p. 407 of Ref. [11]: 


f 


R 2 


OO 


r + ^ n 

n = 2 


(i|)“(c;e3-c; +1 f) 
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(43) 



4.3 Gravitational Force Exerted by an Oblate Spheroid 

An oblate spheroid is an example of an axisymmetric body, and is often used to model the 
Earth. A celestial body’s oblateness is represented by a zonal harmonic of degree 2, and is 
responsible for precessions in a satellite’s orbit plane and argument of perigee. The derived 
Legendre polynomials identified with oblateness are [see Eqs. (6) and (7)] just 

^2,o(sa) = ^( 3s a 2 - 1), ^2,1 (sa) = 3s a (44) 

Substitution from Eqs. (44) and (8) into Eq. (12) yields, for n = 2 and m = 0, 


f 2,0 



3s A e 3 + -(! - 5 s a 



(45) 


This expression for the contribution of oblateness to the gravitational force per unit mass, 
exerted by Earth, is in agreement with the vector equation (8.28) in Ref. [12] and, after a 
little manipulation, with the scalar equations for a x , a y , and a z given on p. 3-1 of Ref. [3]. 
[a x = (-/if /R 2 + f 2 ,o) • ei, ci y = (-/if/R 2 + f 2j0 ) • e 2 , and a z = (-/if / H 2 + f 2;0 ) • e 3 .] 
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Examples show the expressions can be used to reproduce existing results in special cases such as the magnetic field of a tilted dipole, or 
the gravitational field exerted by an oblate spheroid. 
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